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Inhomogeneous plasmas and fluids contain energy stored in inhomogeneity and they naturally tend to relax 
into lower energy states by developing instabilities or by diffusion. But the actual amount of energy in such 
inhomogeneities has remained unknown. In the present work the amount of energy stored in a density gradient 
is calculated for several specific density profiles in a cylindric configuration. This is of practical importance 
for drift wave instability in various plasmas, and in particular in its application in models dealing with the 
heating of solar corona because the instability is accompanied with stochastic heating, so the energy contained 
in inhomogeneity is effectively transformed into heat. It is shown that even for a rather moderate increase 
of the density at the axis in magnetic structures in the corona by a factor 1.5 or 3, the amount of excess 
energy per unit volume stored in such a density gradient becomes several orders of magnitude greater than 
the amount of total energy losses per unit volume (per second) in quiet regions in the corona. Consequently, 
within the life-time of a magnetic structure such energy losses can easily be compensated by the stochastic 
drift wave heating. 

BAGS numbers: 52.35.Kt; 52.50.-b; 96.60.-j; 96.60.P- 


I. INTRODUCTION 


From standard theory it is known that an isothermal 
system with the volume V, consisting of several species 
j and with the total number of particles for each species 
Nj, contains a total amount of internal energy in the 
given volume Uint = Ut = (3kT/ 2) which is just 

the thermal energy of the system. 

This internal energy may be modified in various ways. 
For example, if the species are charged and they include 
electrons and ions, the chaotic thermal motion is affected 
by the Coulomb interaction. The internal energy in this 
case contains such a contribution and it is reduced: 


Uint — Ut 


Utt 



( 1 ) 


Here, d is the mean distance between the charged par¬ 
ticles and ru is the plasma Debye radius. Such a modi¬ 
fication of the internal energy implies also modifications 
of all other thermodynamic functions; i.e., the free en¬ 
ergy, pressure, and entropy are all reduced, and the heat 
capacity increased. One example of this kind, dealing 
with dusty plasmas, may be found in Pandey & Vran- 
jes^^. In most of electron-ion plasmas the Coulomb cor¬ 
rection is small (implying a relatively large amount of 
particles within the Debye sphere), and such plasmas are 
then called ideal, or weakly non-ideal ones. 

The total internal energy (and thus all other thermo¬ 
dynamic functions as well) may change also in the pres¬ 
ence of some external forces. Without such forces a gas 
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or plasma system will tend to relax into an energetically 
most favorable state with minimum energy, and this usu¬ 
ally means it will become isothermal and homogeneous. 
However, these external forces may cause gradients of 
various quantities, like density, temperature etc., and the 
internal energy ([^ and other thermodynamic functions 
will change. Classic studies dealing with the free energy 
and thermodynamics of an inhomogeneous environment 
may be found in works Cahn & Hilliardl^, Hart^^ Rice 
& Chang®, Silver®, McCoy & Davis®, Warner®, and 
Lowett & Bans®. These works mainly contain a general 
theory while for practical purposes some explicit results 
for the free energy are needed, and such an analysis will 
be performed in the present work. 

The presence of such gradients implies an excess of en¬ 
ergy, and the system will tend to get rid of it through 
diffusion of particles until the minimum energy state is 
achieved, or much more efficiently by developing various 
instabilities, including drift instabilities as one example 
of general interest for plasmas. Drift instabilities in an 
inhomogeneous environment represent an efficient way 
of relaxation towards a lower energy state, and they are 
present in almost every realistic plasma configuration. In 
the laboratory plasmas they represent a major issue for 
confinement, though the actual amount of energy stored 
in the plasma inhomogeneity has not been calculated so 
far. For various purposes it may be very useful and in¬ 
structive to know in advance the total amount of this 
extra energy. For example, this may help in modeling 
stochastic heating based on such drift instabilities, in 
particular in the solar corona environment, suggested re¬ 
cently in a series of our works.®^^! 

In the present work the excess energy stored in the 
density gradient is calculated for some relatively simple 
density profiles. The obtained results should be appli- 
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cable to various plasmas both in the laboratory and in 
space, and to fluids in general. 

II. ENERGY IN DENSITY GRADIENT 

We shall assume a cylindric volume V = ttR^L^ inho¬ 
mogeneous due to some external forces, where R is the 
radius and Lz is the axial length. For example, in the 
presence of a magnetic field the radial density gradient, 
or pressure gradient in general, may be counteracted by 
the magnetic pressure. This in principle implies the pres¬ 
ence of the radial inhomogeneity of the magnetic field 
as w ell, w hich can be small in case of a small plasma- 
beta,yet in any case it is supposed to balance the 
assumed pressure (density) gradient. But quite generally, 
the quasi-static equilibrium may be described through 
the balance of forces 

0 = -Vp + Y,Fi- ( 2 ) 

i 

So we may proceed by taking such a density gradient as 
a fact, knowing that it must imply additional forces act¬ 
ing on the plasma, and this means that it is not in the 
minimum energy state. In other words, once such exter¬ 
nal forces are removed, or due to collisions and diffusion, 
the plasma will naturally relax into the lowest possible 
energy state, occupying thus the whole volume V and 
having some constant number density Uh and a constant 
temperature. Our inhomogeneous system thus has some 



r/R 


FIG. 1. Sketch of the two density profiles, linear and Gaus¬ 
sian, in a cylindric volume V = tvR^Lz, in balance with some 
radially acting force F. 

radially dependent density n(r), and we shall calculate 
its internal energy (which includes thermal and poten¬ 
tial energy) taking several possible cases for the density 
shape. 

The theory presented here may be applicable to coro¬ 
nal magnetic structureJi^ [see Fig.[^, and to some atmo¬ 
spheric phenomena like roll clouds and bores. Roll clouds 



FIG. 2. Goronal magnetic structures, picture taken by the 
Transition Region and Goronal Explorer (TRAGE).^^ 



EIG. 3. A roll cloud pictured in Uruguay.^ 



FIG. 4. Multiple roll clouds pictured in Australia.^^ 


are tube-shaped, many kilometers long structures rolled 
up about a horizontal axis, singl^i^ or multiple [see 
Figs. mu]. like so called Morning Glory cloud frequently 
observed in AustralicP^that can be between hundred and 
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thousand kilom eters long and only about one kilometer in 
diameterpUm an d also clouds observed in South Amer¬ 
ica and elsewhere b^ l ^^ l A clear pressure jump of up to 200 
Pa in such rolling solitons has been measuredp^ imply¬ 
ing the excess of energy in them as compared with the 
surrounding fluid. Yet another atmospheric example of 
similar kind are tornados. 

In derivations further in the text, the potential and in¬ 
ternal energies are per species, so the theory is valid for 
ordinary fluids, and in electron-ion plasmas (or multi- 
component in general with a number of species /i) all ex¬ 
pressions should be multiplied by a factor 2 (or by factor 

^l). 

A. Linear density profile with porous boundaries 

We shall assume a nearly ideal plasma with a negligible 
Coulomb interaction energy and having a linear number 
density profile 

n(r) = no — ar. (3) 


with the boundary conditions n(0) = ng, n{R) = Ur this 
yields 

Ut = VKT(^'^+nRy (5) 

In view of Eq. (|^ , the internal energy in the present case 
must be the sum of the thermal energy Ut and some po¬ 
tential energy which system has achieved due to the work 
of the external forces principle, presently we 

are not interested in separate forces Fi, but we shall use 
the fact that their total sum exactly matches the pres¬ 
sure force. So in order to calculate the potential energy 
it may be convenient to use the pressure force instead 
of and express the potential energy through the 

temperature and density only. 

The potential energy Ei(r) of a unit volume at some 
position r, obtained by the act of forces acting on it, is 
the integral over the distance and it can be written as 

''^^Fi{F)dr' = f Vpdr' = KTni{r). (6) 

i 


The number density has the value at the edge of the 
cylinder with the radius i?, and it increases towards the 
axis to a maximum value no, see Fig. the dashed line. 
This means that a = (no — nR)/R. It will be shown later 
that the expected results for the internal energy should 
not be much different from some more realistic cases, like 
the Gaussian profile presented in the same figure. Note 
that in some cases, like in some magnetic configuration in 
the solar plasma, the profile may be just inverted, with a 
decreasing density towards the axis, implying forces act¬ 
ing in the opposite direction. However, there will again 
be a resulting potential energy stored in the inhomogene¬ 
ity. 

If the forces are removed, or due to collisions, in time 
the system will relax to the minimum energy homoge¬ 
neous state with some constant density Uh- In space or 
in solar plasma we may assume that this extra density in 
the gradient from Fig. [^will be absorbed by the infinite 
reservoir of the plasma around with the assumed density 
Ur. So the potential energy will be calculated with re¬ 
spect to riH, and this we shall call the model of porous 
or absorbing boundaries. In such environments this is 
rather justified. However, in bounded plasmas and fluids 
in the laboratory, the relaxed density will clearly have 
some value > Ur. This issue shall be addressed later 


completely exactly in Sec. HI 


The total number of particles in the volume V is 
N = J n(r)dV, where dV = rdrdOdz and integration 
is over the angle 6 S (0,27r), over z € (0, L^), and 
r G (0,R). Eventual variation of density in z and 0 di¬ 
rections is neglected. The result is: 


N = V 


no + 2nj: 


( 4 ) 


The thermal energy for the case of inhomogeneous den¬ 
sity is Ut = f (3KT/2)n(r)dV. For the density § and 


The number density n(r) has the lowest value nR at the 
boundary, and in case of some large configurations in the 
atmosphere or in space, like solar magnetic structures 
this value will be assumed the same outside for r > R 
(porous or absorbing boundary). Inside the structure, 
clearly the excess density above nR will contain addi¬ 
tional potential energy due to action of forces, see Fig. 
Thus the density ni in Eq. (§ satisfies the following 
boundary conditions: ni(0) = ng — nR, and ni{R) = 0. 
This yields ni(r) = (ng — nj,)(l — r/R). 

The total potential energy in the volume V is thus the 
integral of Ei{r) over the whole volume and it reads: 

Ep = VkT -^-. (7) 

The total internal energy is the sum of the two, Uint = 
Ut + Ep, and the Gibbs-Helmholtz thermodynamic free 
energy F can easily be obtained using the well known 
relation 


U =-T^±(l\ 

1 dT\T)y' 

The specific heat and entropy are given by = 
dUint/dT\v, S = —dF/dT\v However, for the drift in¬ 
stabilities only the potential part of the internal energy 
may play the role, and Eq. 0 shows that it is present as 
long as nR ^ ng. 


B. Gaussian and arbitrary flat-top density profile 

For a Gaussian density profile (the full line in Fig. 
n(r) = Cl exp(—ar^) -|- C 2 , where Ci, C 2 are chosen to sat¬ 
isfy the boundary conditions n(0) = ng, n{R) = nR, we 
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have 


C. Application to heating in the solar corona 


( 8 ) 

This is used to calculate the thermal energy. The result 
is 

Ut = ^VKT[a{no - Uji) + n,^], (9) 


a = 1 


1 


1 

1 — exp(—ai?2)' 


For the potential energy the boundary conditions are as 
before ni(0) = uq — n^, ni{R) = 0, i.e., the density Q 
is reduced by rin. The mean potential energy per unit 
volume now becomes 

^ = KTa{no - n^). ( 10 ) 

Comparing this with Eq. Q it is seen that the potential 
energy in the present case may be greater than the value 
obtained for the linear density profile providing that a > 
1/3. Note that a. —)■ 0.5 when —)■ 0, and a —>■ 0 

when aR? —)■ oo, and a ~ 0.42 when aR^ « 1. So the 
differences between the two profiles are not drastic in any 
case. The total internal energy is consequently 

VkT 

U^nt = Ut + Ep = [5a(no - n^) + Sn^]. 


Instead of the Gaussian density we may have an arbitrary 
flat-top profile which can be found in numerous plasmas 
n(r) = Cl exp(—ar^)-|-C 2 , where again ci, C 2 are chosen to 
satisfy n(0) = ng, n{R) = n^. For k = 2 this yields the 
Gaussian case discussed above, and in the limit fc —)■ oo it 
gives a step-profile. Repeating the calculation presented 
above we obtain the potential energy density (101 where 
now 


a = 1 — 


1 

1 — exp(— 


2 1 0 - 2 /'= 
i?2 1 _ exp(—ai?'=) k 


{T[2/k,ar'^])^, 


( 11 ) 


and r[6, a:] is the incomplete Euler gamma function. 

In fact, the analysis can easily be generalized to any 
density profile in the cylindric system, bearing in mind 
the procedure described above. The result for the poten¬ 
tial energy is 

2kT 

= ^~j^ J [”(^) “ ( 12 ) 

and the thermal energy is Ut = SttkTLz n{r)rdr = 
{VZkT/R?) n(r)rdr. 


The solar corona is hot and it continuously looses its 
energy. Without an efficient generator of energy, the 
amount expressed through Eq. Q would be lost in less 
than a day and the corona would become completely cool. 
However, it remains hot and with the temperature ex¬ 
ceeding a million degrees. 

The presence of the density gradient implies drift in¬ 
stabilities and those may lead to stochastic heating if 
certain conditions are satisfied.^^^^^^ We shall make some 
estimate and application of the present results to the drift 
wave heating paradimn in the solar atmosphere, which we 
put forward recently.!®!^ 

For the quiet regions in the solar corona, presently 
accepted required heating rate,l^ with the temperature 
T = 1.8 • 10® K, is in the range Fg = Eg/s = 2 • 10“^ 
J/(m®s) to Fg = 4 • 10“® J/(m®s) for the two number 
densities Uq = 0.2 • lO'^® m”® and Uq = 10®® m”®, respec¬ 
tively. Note that for such parameters the corrections due 
to Coulomb interaction in Eq. Q are completely negli¬ 
gible. Such a heating rate can easily be satisfied in the 
solar magnetic structures following the models based on 
the drift wave stochastic heating.lSHlllln these studies var¬ 
ious heating rate values are obtained, dependent on the 
plasma parameters. As example, in magnetic structures 
with the characteristic inhomogeneity scale-length L„ of 
around several hundred kilometers, the wave frequency is 
of the order of 0.1 Hz and the energy release rate F due 
to the stochastic heating is around the value Fg given 
above, see more in Vranjes & Poedts^^. 

We may now take some specific density profile and cal¬ 
culate the potential energy density Ei = Ep/V for the 
heating rate F ~ Fg in order to see at least roughly for 
how long such a density profile might sustain the plasma 
at the given temperature, i.e., for how long it may com¬ 
pensate for the energy losses. The result is presented in 
Table |T] and Table [IT] for the temperature and two start¬ 
ing densities for quiet regions Un = Uq given above, and 
allowing for several possible density values ng in the cen¬ 
ter of the loop. From Table |T] it is seen that even in the 
case when the boundary value increases only by a 
factor 1.5 at the axis, the mean potential energy density 
becomes Ei = Ep/V = 8.3 • 10“"® J/m®, which is several 
orders of magnitude greater than the energy lost per unit 
volume in one second. This is even more so for the values 
in Table [Hj But in reality, the density gradient changes 
(i.e., the density profile flattens) in time and so does the 
realistic energy release rate; observations show that the 
mean life time of the loops is in the range of half an hour 
to a few hours. So to make some estimate we may set 
the loop life-time equal one hour, and to have a realis¬ 
tic margin and to be sure about a sustainable heating 
it would be necessary to have the ratio Ei/Eg consider¬ 
ably greater than 3600. The numbers given in Tables [ij 
ITTl show that this is indeed the case for almost all values 
of ng. Observe also that the potential energy becomes a 
considerable part of the total internal energy already for 
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TABLE I. Parameters for quiet regions following Narain & 
Ulmschenider^ with Ur = 0.2 • 10^® and for linear den¬ 
sity profile given by Eq. with several values no at the 
center. 


ml Ur 

El [J/m®] 
Uth/V [J/m®] 
Ep/Uth 

Ei/Eo 

1.5 

0.0008 

0.0087 

0.095 

4142 

3 

0.0033 

0.012 

0.27 

16568 

5 

0.0066 

0.017 

0.38 

33137 

10 

0.015 

0.03 

0.5 

74558 

TABLE 11. Parameters for 
Ulmscheniderl^ with Ur = 
profile 1^. 

quiet regions following Narain & 
10^® and for linear density 

mjnR 

1.5 

3 

5 

10 

El [J/m®] 

0.0041 

0.0165 

0.033 

0.074 

Uth/V [J/m®] 

0.043 

0.062 

0.087 

0.15 

Ep/Uth 

0.095 

0.27 

0.38 

0.5 

E\/Eo 

20710 

82842 

165684 

372789 


nRlriQ > 3. As stressed before, the numbers presented 
here are per species. 

A similar analysis can be performed for the coronal 
holes where Fq = 8 • 10“® J/(m^s), showing that the 
amount of energy in the density gradient is more than 
enough for a sustainable heating. 

In active regions, following Narain & Ulmscheniderl^ 
the temperature and the number density vary from T = 
10^ K, Ua = 0.5 • m-3 to T = 2.5 • 10® K, «„ = 5 • 

m“^ = Ur, and the energy losses are in the range 
Fo = 7 • lO-'^ J/(m3s) to Fo = 3 • lO"'* J/(m3s). For 
the first set of data a sustainable heating F ~ Fq appears 
possible for a few second only, provided that no = lOn^. 
However, for the second set T = 2.5 • 10® K, n^ = 5 • 10^® 
m“® the situation is a bit different, and the result is pre¬ 
sented in Table EH These data indicate that a rather 
strong density gradient is needed to sustain the energy 
losses within the life-time of a magnetic structure in ac¬ 
tive regions. 

The analysis presented above can be repeated for the 


Gaussian density profile M. Without going into detail 
we may take the case al^ « 1, and in view of Eqs. (H 


101 it may be shown that the above given values for the 


heating time should be multiplied by a factor 1.26. So 


TABLE III. Parameters for active regions following Narain & 
Ulmschenider^ with T = 2.5 • 10® K, ria = 5 • m“®, for 

linear density profile §, and with several values no at the 
center. 


m/nR 

1.5 

3 

5 

10 

15 

El = Ep/V [J/m®] 

0.03 

0.115 

0.23 

0.52 

0.18 

Uth/V [J/m®] 

0.3 

0.43 

0.6 

1.04 

1.47 

Ep/Uth 

0.095 

0.27 

0.38 

0.5 

0.55 

Ei/Eo 

96 

383 

767 

1726 

2685 


the heating in this case in quiet regions and coronal holes 
is even more certain and this holds for all cases with 
aR^ < 1 . 

The calculated potential energy stored in the density 
gradient and the consequent stochastic heating by the 
drift wave are thus more than enough to compensate 
for losses in quiet Sun regions and in coronal holes, and 
the predicted sustainable stochastic heatinJ^HSl by the 
drift wave looks like a rather realistic scenario. In active 
regions it may sustain heating only provided relatively 
strong density gradients as Table m suggests, so most 
likely some additional and more energetic processes are 
in action in such regions. 


III. NON-ABSORBING BOUNDARIES 


In the laboratory environment the initial density pro¬ 
file with the number density Ur at the edge of the cylinder 
will normally relax to a homogeneous minimum available 
energy state with the density line n = nt, see Fig.[^ We 
shall use this simple linear density for the present analy¬ 
sis although obviously it can easily be replaced with the 
Gaussian one as shown in Sec. |IT] Quantitatively, the re¬ 
sults will be similar. Clearly there must be Ur < nu < no 
as depicted in Fig.[^ The same situation may sometimes 
happen also in coronal magnetic structures with excep¬ 
tionally strong boundary magnetic field when diffusion 
across the boundary itself is much weaker as compared 
with the diffusion inside the structure and when there are 
drift instabilities taking place inside the structure. In this 
case the density nu would correspond to the density of the 
bulk plasma around the structure. A more realistic sit¬ 
uation in coronal structures is presented in Fig. where 
we have the external density n^, and the given density 
inside the structure, which is a more physical variant of 
the rough profile from Fig. 

Note that the assumed actual linear density profile in 
Fig. [5] is in fact partly inverted with respect to nh, and 
both domains n(r) > Uh and n(r) < have some poten¬ 
tial energy with respect to the relaxed state with n = Uh- 

The density can be calculated in the following man¬ 
ner. The total number of particles in the volume V is 
again given by Eq. Q. In the relaxed state, both the 
total number of particles and volume remain the same, 
so N = UhV and this directly yields 

no , 2nR 

The intersection of the lines n = n{r) and n = nh 
at Tc = 2i?/3 and this remains so for any no and Ur. 
At t —)■ 00 the internal energy will reduce to the thermal 
energy In any other moment it will be the sum of the 
internal energy (j^ and potential energy Ep with respect 
to the density level nh- The potential energy is the sum 
of the two integrals Sp = Epi -|- Sp 2 , 


Epi = kT 


[n(r)-nh]dV, 


Sp 2 = kT f [nh-n{r)]dV. 


0 
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FIG. 5. Sketch of the density profile in a cylindric volume 
V = tvR^Lz for a non-absorbing boundary. Without forces 
the system relaxes to a homogeneous state with n = nu where 
Ur < Uh < no. 


FIG. 6. Sketch of the two density profiles, linear and Gaus¬ 
sian, in a cylindric volume V = tvR^Lz, in balance with some 
radially acting force F. Without forces the system relaxes to 
a homogeneous state with n = Uh- 


The result turns out to be 

g 

Spi = T.p2 = —VkT^uo - Ur), 

so that 

T,p=^VKT{no-nR). (13) 

This may be compared with the potential energy Q from 
the case of absorbing boundaries. The ratio of the two 
is Ep/T,p « 1.7. Such a result could have been foreseen 
because the reference level with respect to which the en¬ 
ergy is calculated in the present case is higher, so the 
potential energy is naturally smaller. 

It should be stressed that this is a completely general 
and exact analysis and it can be performed in the same 
manner for any other analytically given density profile. 
For example, the profile in Fig. [^is in fact a function of 
the kind a + exp(—C 2 r^)Jo(c 3 r) (where Jg is the Bessel 
function), but such derivations will not be repeated here. 

IV. SUMMARY AND CONCLUSIONS 

Inhomogeneity of plasmas implies a free energy for de¬ 
velopment of drift instabilities as the most efficient way 
of relaxing into a more favorable lower energy state. Such 
instabilities may be accompanied with stochastic heating 
provided large enough amplitudes of the perturbations, 
and this phenomenon has been observed in the past in 
specially designed laboratory experiments.l^^l^^ 

The same idea has been used recentlj^tHl as a new 
paradigm for the heating of the solar corona. Analytical 
calculations in these references show that such a stochas¬ 
tic heating can easily satisfy numerous heating require¬ 
ments like the necessary heating rate, better heating of 


heavier particles, stronger heating in perpendicular direc¬ 
tion with respect to the magnetic field vector, etc. The 
model is based on a quasi-static inhomogeneous back¬ 
ground plasma which is supposed to contain enough en¬ 
ergy in order for the mechanism to work. This assump¬ 
tion is checked in the present study, and the energy con¬ 
tained in the density gradient is calculated. The results 
obtained in the work suggest that accidentally and exter¬ 
nally created magnetic structures indeed contain enough 
energy for a sustainable heating at least in the solar quiet 
regions and in coronal holes. 

The inhomogeneity is an intrinsic feature of the mag¬ 
netic structures in the solar atmosphere. Movement and 
restructuring is a continuous process in coronal magnetic 
structures and such phenomena are accompanied with 
the motion of plasma due to its frozen-in properties. The 
life-time of these structures is of the order of an hour, so 
that indeed they can be treated as quasi-static for rel¬ 
atively fast drift instabilities, and they are maintained 
by processes deep below the corona itself. This implies 
that such a drift-wave-favorable environment is being cre¬ 
ated continuously; structures appear and disappear all 
the time, and heating is taking place in each of them in¬ 
dependently. As discussed in Vranjes & Poedts^^Eni^ the 
drift waves and heating are produced directly in these 
magnetic structures and driven by the plasma inhomo¬ 
geneity. However, the complete mechanism is totally 
magnetic by nature: such inhomogeneities are caused by 
the dynamics and restructuring of the magnetic field, and 
these phenomena on the other hand are driven by pro¬ 
cesses far from the corona, i.e., inside the Sun. In the 
present work we were focused only on the heating rate 
close to the energy loss rate, to check if it can be sustain¬ 
able or not. However various heating rates are possible 
dependent on the drift wave frequency and growth-rate 
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(which on the other hand are dependent on the inho¬ 
mogeneity scale-length), including energy releases in the 
range of nano-flares^. 

The analysis presented in the work is for a simple cylin- 
dric geometry, but it can easily be generalized to other ge¬ 
ometries applicable to the laboratory plasma like a toka- 
mak, or to loop structures in the solar atmosphere. 
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